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After a quench, observables in an integrable system may not relax to the standard thermal values, 
but can relax to the ones predicted by the Generalized Gibbs Ensemble (GGE) [M. Rigol et at, 
PRL 98, 050405 (2007)]. The GGE has been shown to accurately describe observables in various 
one-dimensional integrable systems, but the origin of its success is not fully understood. Here 
we introduce a microcanonical version of the GGE and provide a justification of the GGE based 
on a generalized interpretation of the eigenstate thermalization hypothesis, which was previously 
introduced to explain thermalization of nonintegrable systems. We study relaxation after a quench 
of one-dimensional hard-core bosons in an optical lattice. Exact numerical calculations for up to 10 
particles on 50 lattice sites (sa 10 10 eigenstates) validate our approach. 

PACS numbers: 02.30.Ik,03.75.Kk,05.30. Jp,67.85.Hj 
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Once only of theoretical interest, integrable models of 
one-dimensional (ID) quantum many-body systems can 
now be realized with ultracold atoms 0. The possibility 
of controlling the effective dimensionality and the degree 
of isolation have allowed access to the quasi- ID regime 
and to the long coherence times necessary to realize in- 
tegrable models. Additionally, advances in the cooling 
and trapping of atoms have led to increased interest in 
dynamics following quantum quenches, where a many- 
body system in equilibrium is exposed to rapid changes 
in the confining potential or interparticle interactions. 

In general, in integrable quantum systems that are far 
from equilibrium, observables cannot relax to the usual 
thermal state predictions because they are constrained 
by the non-trivial set of conserved quantities that make 
the system integrable 0] . Relaxation to non-thermal val- 
ues were recently observed in a cold- atom system close to 
integrability Q . At integrability, it is natural to describe 
the observables after relaxation by an updated statisti- 
cal mechanical ensemble: the generalized Gibbs ensem- 
ble (GGE) 0, which is constructed by maximizing the 
entropy subject to the integrability constraints |5|. In 
recent studies of integrable systems [1, 0, 0|, the GGE 
has been found to accurately describe various observables 
after relaxation, but a microscopic understanding of its 
origin and applicability remains elusive. In particular, an 
important question remains: how is it that expectation 
values after relaxation can be described by an ensem- 
ble with exponentially fewer parameters than the size of 
the Hilbert space? The full dynamics are determined by 
as many parameters as the size of the latter. At a mi- 
croscopic level, thermalization for non-integrable systems 
can be understood in terms of the eigenstate thermal- 
ization hypothesis (ETH) [1,11], which, however, breaks 
down as one approaches integrability [lOj. 

This paper is devoted to the study of how generalized 
thermalization, in the sense of relaxation to the predic- 
tions of the GGE, takes place in integrable systems. An- 



swering this question is important not merely because of 
its relevance to the foundations of statistical mechanics in 
integrable systems, but also because it has become nec- 
essary to understand recent experiments with ultracold 
gases in quasi- ID geometries. For integrable systems, 
we compare the predictions of quantum mechanics with 
those of various statistical ensembles. In particular, we 
introduce a microcanonical version of the GGE, which 
we use to show that relaxation to the GGE can be un- 
derstood in terms of a generalized view of the ETH. 

We study the dynamics following an instantaneous 
quench of ID hard-core bosons on a lattice, which is fully 
integrable. The Hamiltonian is given by 

L-l L 

H = -jJ2 (Hh+i + H.c.) + V(t) ^(i- L/2) 2 h l (1) 



where J is the hopping parameter; V(r) gives the cur- 
vature of an additional parabolic trapping potential for 
atoms on a lattice with lattice constant a; b\ (pi) is the 
hard-core bosonic creation (annihilation) operator; and 
hi = b\bi is the number operator. In addition to the stan- 
dard commutation relations for bosons, hard-core bosons 
satisfy the constraint b\ 2 = tij = 0, which forbids mul- 
tiple occupancy of the lattice sites. This Hamiltonian 
can be mapped onto non-interacting fermions through 



the Jordan- Wigner transformation and the many- 
body eigenstates can be constructed as Slater determi- 
nants of the single-particle fermionic eigenstates [lj . 

We will focus on the behavior of the momentum dis- 
tribution function, (%) = Ei,m e ' ik{l ~ m) ('>PP m bi\^}/L, 
for system sizes ranging from N = 5 bosons on L = 25 
lattice sites to N = 10 bosons on L = 50 lattice sites 
(s=s 10 10 eigenstates). Initially, we prepare the system 
in the ground state |^o) of a ID lattice with hard- wall 
boundary conditions and an additional harmonic poten- 
tial, with trapping strength V = Vq. At time r = 0, the 
harmonic trap is turned off, V(t > 0) = 0, and the state 
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FIG. 1. (a),(b) Momentum distribution of the initial state 
(init), diagonal (DE), generalized microcanonical (GME), 
generalized Gibbs (GGE), and the microcanonical (ME) en- 
sembles. (c),(d) Relative difference of the GME, GGE and 
ME from the DE. (e),(f) Conserved quantities, (In), in the 
quenched state (identical to the DE), GME and ME. (J„) 
are ordered in descending occupations in the quenched state. 
L = 50, TV = 10, 5me = 0.05J, <5 G me = 0.8. (a),(c),(e) 
e = 0.72J, V = 0.029J. (b),(d),(f) e = 1.52J, Vo = 0.125J. 

\ip(r)) evolves under the influence of the final Hamilto- 
nian. Hereafter, we refer to this state as it is immediately 
after the quench as the "quenched state". Its time evo- 
lution is given by \ip{r)) — ^2 a c a e~ lEaT l h \a) , where |a) 
are the energy eigenstates of the final Hamiltonian with 
energies E a , and c Q = (a\ipo) are the overlaps between 
the eigenstates of the final Hamiltonian and the quenched 
state. After relaxation, assuming the degeneracies in en- 
ergy levels are irrelevant, the expectation value of an ob- 
servable is expected to be given by the so called diagonal 
ensemble (DE) @, 11 E3 
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We have checked numerically that, despite the integra- 
bility of our model, n k relaxes to the DE prediction, with 
small fluctuations around this result fl3j |. 

Figure 1 shows the momentum distributions, rife, be- 
fore and after the quench, for two different initial trap 
strengths, which correspond to different energies per par- 
ticle, e, after the quench. These results are compared 
with those of various ensembles of statistical mechanics. 
The microcanonical ensemble (ME) is one in which all 
eigenstates in the relevant energy window have identi- 
cal weights. Within the microcanonical ensemble, the 
expectation value of a generic observable A is (A)me = 

N f ;i ME E a j e -e Q |<5 ME H^l a >' where Sme is sma11 ' but 

still much greater than the mean many-body level spac- 
ing. NsJmb i s the number of eigenstates in the energy 



window |e — e a \ < Sue- We have checked that the re- 
sults reported here are nearly independent of the specific 
value of <5me- The GGE is a grand-canonical ensemble 
that maximizes the entropy subject to the constraints 
associated with non-trivial conserved quantities of the 
quenched state. The density matrix takes the form Q 



PGGE = Z G 1 e Z G = Tr 



(2) 



where {/„}, n = 1,...,L, are the conserved quanti- 
ties. In our systems, these correspond to the occupa- 
tion of the single-particle eigenstates of the underlying 
nonintcracting fermions to which hard-core bosons are 
mapped, and {A„} are Lagrange multipliers fixed by the 
initial conditions, A„ = ln[(l - (ip \I n \ip }) / {if} \I n \ipo}} 
[3| . Observables within this ensemble are then computed 

as (A)qge = Tr Ap GGE following Ref. 

As a step towards understanding the GGE as well as 
developing a more accurate description of isolated inte- 
grable systems after relaxation, we introduce a micro- 
canonical version of the GGE, which we call the gen- 
eralized microcanonical ensemble (GME). Like the ME, 
where states within a small energy window contribute 
with equal weight, within the GME we assign equal 
weight to all eigenstates whose values of the conserved 
quantities are close to the desired values. The expec- 
tation value of a generic observable within the gener- 
alized microcanonical ensemble is given by (^4)gme = 
AT/i x y~!„ x (oi\A\a), where Y\ A is a 

sum over eigenstates that are within the GME window 
and ■A/{/„} ) 5 GME is the number of states within that win- 
dow and 8 a is a measure of the distance of eigenstate a 
from the target distribution. 

In order to construct the GME, we include eigenstates 
of the Hamiltonian with a similar distribution of con- 
served quantities which once averaged reproduce the val- 
ues of the conserved quantities in the quenched state. 
This approach is characterized by three ingredients: (i) 
The ordered distribution (from largest to smallest) of the 
conserved quantities in the DE, (I„)de = J2a I 
[as in Figs. 1(e) and 1(f)], (ii) a target distribution of 
the nonzero expectation values of the conserved quanti- 
ties {I*, =1}, where the values of n* (i — 1, . . . , TV) are 
chose to describe the distribution /„ in a coarse grained 
sense [l5| , and (iii) for each individual many-body eigen- 
state, the distance from the target state, <5 Q , which we 



define as 8 a = 



1/2 



Here , 



(i = 1,...,TV) are the single-particle states occupied 
in eigenstate a, and /„* are the interpolated values of 
(^ti)de, evaluated at n*. The definition of 5 a is not 
unique and several variants that do not change our con- 
clusions were also considered [131 ] . 

To better visualize the differences between the re- 
sults of the various ensembles in Figs. 1(a) and 1(b), 
we have plotted A(n k ) sta t = ((%)de - (%)stat)/ («)c)de, 
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where "stat" stands for ME, GGE, or GME in Figs. 1(c) 
and 1(d). For weaker initial confinements (smaller e - 
Fig. 1(c)), the GME is practically indistinguishable from 
the diagonal distribution. Both the GME and the GGE 
accurately capture the tails of rik, while the thermal en- 
semble does not. For tighter initial traps (greater e - 
Fig. 1(d)) all four ensembles are very similar (note the 
scale), suggesting that rik in the final steady state is in- 
distinguishable from that of the thermal state. 

The close agreement between DE and ME results in 
Fig. 1(b) raises the question: how can an integrable sys- 
tem thermalize, given the constraints imposed by the 
complete set of conserved quantities? We conjecture that 
if the values of the conserved quantities in the quenched 
state are similar to those of the ME, then the latter will 
accurately describe observables after relaxation. This 
may occur for a variety of quenches. 

In Figs. 1(e) and 1(f), we plot the values of the con- 
served quantities in the quenched state and compare 
them with the expectation values of those quantities in 
different statistical ensembles. (By definition, the dis- 
tribution of conserved quantities in the DE and GGE 
are identical to that of the quenched state.) Figure 1(e) 
shows that the microcanonical values of the conserved 
quantities are clearly different from the values in the 
quenched state, while in Fig. 1(f) they are very similar. 
This supports our conjecture above, and demonstrates 
that thcrmalization can occur in integrable systems for 
special initial conditions. Additionally, the GME repro- 
duces the correct distribution of the conserved quantities 
supporting the validity of our method for generating it. 

To quantify the above observations, and to under- 
stand what happens in the thermodynamic limit, we 
have studied the difference between the predictions of 
the DE and the statistical ensembles for different system 
sizes. We compute the integrated relative differences, 

(Afjfc)stat = 2feK%)DE - («fc)stat|/Z)fc(n*i)DE, where 
again "stat" stands for ME, GGE, or GME. 

In Fig. 2(a), we plot (Aii^Jme as a function of the 
final energy per particle, e, for different lattice sizes, 
L. To perform finite-size scaling, e and the filling fac- 
tor (y — N/L — 0.2) are held constant as L changes. 
Figure 2(a) shows that for e < 1.3J the difference be- 
tween the rik in the DE and the ME increases with 
increasing L, indicating that the difference persists in 
the thermodynamic limit. For e > 1.3 J, the oppo- 
site behavior is observed. From our previous discus- 
sion, one expects that (Ati^me should closely follow 
the behavior of the integrated differences between the 
conserved quantities in the quenched state and the ME, 

(A/„) M E = S n |(I«.}de - (^)me|/IZ„(-^)de- This is 
seen by comparing Figs. 2(a) and 2(c), which leads us 
to conclude that rik need not relax to the standard ther- 
mal prediction, except when (A/„)me becomes negligi- 
ble. Qualitatively similar results were obtained in the 
canonical ensemble 1131 . 
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quenched state. <5me = 0.05J. (b) (An fc )GME vs e, 5gme = 
0.8. (c) Integrated difference between the conserved quan- 
tities in the quenched state and the ME, (AJ„)me- (d) 
(Anfc)GGE vs e. Inset: (An fc )GGE vs L -0 ' 73 for e = 1.07J, 
where a fit to (Arife)GGE = zL' 1 gives 7 = 0.73 ± 0.02. 



On the other hand, in Fig. 2(b) one can see that the dif- 
ferences between rik in the diagonal and generalized mi- 
crocanonical ensembles are very small and decrease with 
increasing system size, so that the former successfully 
describes this observable after relaxation. In the case of 
the GGE [Fig. 2(d)], (Anfc)cGE is in general larger than 
(An fc ) G ME, which is to be expected since the GGE is a 
grand-canonical ensemble. As the system size increases 
(An fe ) GGE -> as L"T, where 7 « 0.73 [inset of Fig. 2(d)] 
and slightly depends on the energy [lj|. 

The question that remains to be answered is why the 
generalized Gibbs and the generalized microcanonical en- 
semble are able to describe the rik after relaxation, i.e., 
why (n^QGE = ("-Ogme = ("-Ode = J2 a \c a \ 2 {a\h k \a) . 
Note that whereas (nfc)GGE and (b^gme are entirely de- 
termined by the L independent values of the conserved 
quantities in the quenched state, (nfc)DE is determined 
by the exponentially larger (^j.) values of the coefficients 
c a . 

To address this question, we perform a spectral de- 
composition of (fifc)DE and (rik) gme- Figure 3 displays 
a coarse grained view of the weight which eigenstates 
with a given zero momentum occupancy (hk=o)a — 
(a|nfc=o|a) contribute to the DE [Fig. 3(a)] and the GME 
[Fig. 3(b)]. The correlation between the results in both 
figures is apparent. However, it is not clear why the de- 
tails contained in the overlaps c a are completely washed 
out so that the DE and the GME results coincide, while 
they are different from those in the ME. In the inset of 
Fig. 3(a), we plot a histogram of the values of Hfc = o f° r 
the DE, GME and the ME. Clearly the histograms for the 
DE and GME have a similar mean but different widths, 
while the ME has a different mean and width (l3| . 

Ultimately, one is interested in what happens in the 
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FIG. 3. Density plot of the coarse-grained weights with 
which eigenstates contribute to (a) the DE (sum over diagonal 
weights, |c a | 2 ) and (b) the GME (fractional number of states 
= number of states/total number of states) as a function of 
eigenstate energy and {hk=o)a- The sums are performed over 
window of width Sn k = 0.0067 , and 5e = 0.035J. The hori- 
zontal and vertical dotted lines are the expectation values of 
nfe = o and e in each ensemble. L — 45, N = 9, Vq = 0.036J, 
e — 0.72J, 5gme = 0.85. Inset in (a): Histogram of DE 
weights (green), fractional number of GME states (blue) and 
fractional number of ME states (black) summed over all ener- 
gies. Bin width, Snt = 0.0067. Vertical lines give the mean, 
{hk=o) within each ensemble. Inset in (b): Fluctuations 
of (n fea=0 >(«) and {n ka =2ir/s){*) within the DE (green) and 
GME (blue) as a function of inverse system size, e = 0.72J. 



thermodynamic limit, 
of the distribution of 



For each k, we define the width 
(nk)a f° r each ensemble as Ofc = 
(fife) 2 . The inset of Fig. 3(b), shows Cfe within 
the DE and the GME versus The scaling is de- 

picted for two k values and clearly shows that the widths 
of both distributions vanish in the thermodynamic limit. 
This demonstrates that the overwhelming majority of the 
states selected by the DE as well as by the GME, which 
have similar values of the conserved quantities, have iden- 
tical expectation values of ■ This is why details of the 
distribution of c a no longer matter as L increases. We 
note that with increasing L, the number of eigenstates 
contained in the generalized microcanonical window in- 
creases exponentially, however, the ratio of the number 
of states in the GME and the ME vanishes 



13] 



The findings above provide a generalization of the ETH 
introduced previously to understand thermalization in 
nonintegrable systems [1, Q . The ETH states that the 
expectation values of few-body observables in generic sys- 
tems do not fluctuate between eigenstates that are close 
in energy. Thus all eigenstates within a microcanoni- 
cal window have essentially the same expectation values 
of the observables, and one can say that thermalization 
occurs at the level of eigenstates. As seen in Fig. 3, 



(fife) a exhibits large eigenstate-to-eigenstate fluctuations 
in our integrable system, showing that ETH is invalid. 
However, by selecting eigenstates with similar conserved 
quantities, ETH is restored, although in a weaker sense: 
the overwhelming majority of eigenstates with similar 
conserved quantities have similar values of rife. These 
results pave the way to a unified understanding of ther- 
malization in generic (nonintegrable systems) and its gen- 
eralization in integrable systems. This opens many new 
questions, such as whether the concepts of typicality [3] 
and thermodynamics 17J can be generalized to isolated 
integrable systems. 
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Time Dynamics. In order to check that the DE ac- 
curately describes observables after relaxation, we com- 
pare the time dynamics of the central momentum peak, 
(fT-fc=o( r )) = (V , ( T )l"-fe=o|^'(''")), with the expectation 
value of hk=o in the diagonal ensemble. Figure 1(a) 
shows that (hk=o(T~)) relaxes to the diagonal prediction, 
with small fluctuations around this result, indicating that 
the diagonal ensemble correctly predicts the values of ob- 
servables after relaxation. Additionally, the expectation 
value of hk=a in the GME agrees with the DE, while the 
ME does not. 

In order to understand the relation between the 
DE prediction and the actual time average over a fi- 
nite time interval, nj: = T Ji ^ dr {h^(r)) , we study 
the integrated difference (Anfc) T = J^^K^a^de — 
Sfe("fe)rJE, and the time fluctuations ay = 
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Sfc \/ n fe — n k 2 as a function of system size. (Ant) T is 
depicted in Fig. 1(b), where we also show the results 
of a fit to (An fc ) T = zL~i with 7 = 0.98 ± 0.25 (z is 
a generic multiplicative constant). The integrated dif- 
ference, (An/s) T , clearly decreases with increasing sys- 
tem size and is several orders of magnitude smaller than 
similar comparisons with the statistical ensembles, con- 
firming that degeneracies are irrelevant and the diagonal 
distribution accurately represents the long-time average. 
Figure 1(c) depicts the time fluctuations of the momen- 
tum distribution, ay , where the exponent was determined 
from fitting ay = zL 1 . This plots makes evident that 
the fluctuations also decrease with increasing system size 
(likely oc 1/VZ) and are expected to vanish in the ther- 
modynamic limit. 

Constructing the generalized microcanonical ensemble. 
We expand upon our method for generating the gen- 
eralized microcanonical ensemble with a series of plots 
in Fig. 2. Constructing the generalized microcanonical 
ensemble for a system of ID hard-core bosons presents 
some unique challenges. In particular, the set of con- 
served quantities are the single-particle occupations of 
the underlying fermions used to build the many body 
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FIG. 1. (a) Time evolution of nk=o- Horizontal lines represent 
<n fc ) D E, (nfe)GME, (nk}ME. L = 50, N = 10, e = 0.72J, 
5gme = 0.8, 5me = 0.05J. Main Panel: r < 100. Inset: 
100 < r < 1000. (b) Integrated difference of the diagonal and 
time-averaged momentum distribution, (Anfc) T , versus L -0,98 , 
where a fit to o T — zL -1 gives 7 = 0.98 ± 0.25. (c) Mean 
fluctuations of the momentum distribution, a T , versus L -0 ' 4 , 
where a fit to (Ara t ) T = zh' 1 gives 7 = 0.40 ± 0.12. (b),(c) 
Ti = 100J, T2 = 1000J, e = 0.72J. 2 is a generic multiplicative 
constant. 



FIG. 2. (a) Distribution of conserved quantities in the diago- 
nal ensemble, (J„)de, ordered by energy of the single-particle 
fermionic eigenstates. (b) (/„)de sorted in descending order. 
The target distribution of n* are labeled as well as the corre- 
sponding J„« for i = 1,2. (c) Target distribution of conserved 
quantities, 7^, and distribution of eigenstate that minimized 
the distance 8 a used to generate the GME, I n , a . (d) (J n ) in 
the DE and GME using method described in text as well as 
the standard metric (GME-std). L = 20, N = 4,e = 0.72J. 
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eigenstates. Thus in any given eigenstate, each conserved 
quantity is either or 1 corresponding to whether or not 
that fermionic state is occupied, such as I n , a in Fig. 2(c). 
On the other hand, the distribution in the diagonal en- 
semble is a continuous variable, < (I„)de < 1- For 
each individual eigenstate, we must include or exclude it 
from the GME based on how close the discrete distribu- 
tion in the eigenstate is to the continuous distribution 
of the quenched state. One possible method is a simple 
extension of the microcanonical metric 



£cc-/»,«) 3 /E 



n=l 



1/2 



We found this was not guaranteed to accurately de- 
scribe the distribution of conserved quantities, particu- 
larly for low energies for the system sizes studied (see Fig. 
2(d)). Instead of using the standard metric, we perform 
a weighted least-squares fit to a target distribution which 
is described as follows. 

The unsorted conserved quantities of the quenched 
state (J„)de, ordered by energy of the single-particle 
fermionic eigenstates, are plotted in Fig. 2(a). As a 
first step, the distribution of (/„)de is sorted in descend- 
ing order as shown in Fig 2(b). Next the values of n*, 
which are not limited to integer values, are chosen so that 
J„ g I(x) = 0.5, where I(x) — (I„)de for x in the interval 
(n — 0.5, n + 0.5]. Subsequently the n* are determined so 
Qi I{x) = 1 as depicted in Fig. 2(b). This set of {n*} 

determine the target distribution, labeled /* in Fig. 2(c). 
Each n* is assigned a weight, /„* , which is the interpo- 
lated value of (J„)de at n*. The distance between the 
distribution of conserved quantities in eigenstate a and 
the distribution of the quenched state is then defined as 



8 a = 



1 N 



1/2 



(1) 



where the the single particle eigenstates that are 

occupied in the many-body eigenstate a. 

Our method of constructing the GME is not unique. 
We tested various different approaches, including an un- 
weighted sum in Eq. ([1]), etc, which gave similar re- 
sults. We then employed the method which best recon- 
structed the distribution of conserved quantities over the 
full range of parameters studied. Note that for clarity 
Fig. 2 displays data for the case of N — 4 particles on 
L = 20 sites. There are only 4845 total eigenstates and 
174 eigenstates used to construct the GME for the data 
shown. Even for such a small system, the distribution of 
conserved quantities within the GME is quite good. 

Lagrange multipliers and additivity. In Fig. 3, we 
plot the expectation values of the conserved quantities 
in the diagonal ensemble along with corresponding La- 
grange multipliers. The Lagrange multipliers are given 



by A n = In ( 1 J J "p° E J ■ As can be seen, the Lagrange 
multipliers vary smoothly with the value of I n and with 
the index n after the conserved quantities have been or- 
dered in descending order. Additionally, the distribution 
of conserved quantities is very similar for two different 
lattice sizes with the same final energy-per-particle, when 
the index is normalized by the total number of conserved 
quantities. This is an important property of the Lagrange 
multipliers, which shows that even though the conserved 
quantities in these integrable systems are not additive in 
a strict sense, they can still be understood to be additive 
in a coarse grained sense, because the values of X n are a 
smooth function of /„ . 
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FIG. 3. Expectation values of the conserved quantities in the 
diagonal ensemble and corresponding Lagrange multipliers for 
N = 5, L = 25 and N = 10, L = 50. e = 0.72i 

Canonical Ensemble. In addition to the results re- 
ported for the microcanonical distribution, we study 
the momentum distribution in the canonical ensemble 
(CE). Within the canonical ensemble, observables are 



calculated as (A)ce 



Ae-W 



where Zce 



Tr 



and /3 is the inverse temperature. The tem- 
perature is calculated numerically so that (E)ce — 
(ipo\H(r > 0)|^o)- In Fig- 4(a) we plot the integrated 
difference between the momentum distribution in the DE 
and CE, (Anfc)cE, for different lattice sizes as a func- 
tion of the energy per particle of the quenched state. 
The comparison of the canonical and diagonal momen- 
tum distributions are similar to the comparison between 
the microcanonical and diagonal distributions, although 
there is some discrepancy between the microcanonical 
and canonical distributions due to finite size effects. In 
particular, the upturn in (An/^ME at e — 1.4J is not 
present in (Att^ce- 

In general, the microcanonical distribution provides us 
with a better description of the system because it is iso- 
lated, although it breaks down for very small systems 
because of poor statistics due to an insufficient num- 
ber of states in the relevant energy window. Given 
the expectation that (hk)cE and (n^ME are equal in 
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of the exponent is 7 = 0.731 ± 0.007. 



0.01 0.02 

I -1-0 



FIG. 4. (a) Integrated difference between the diagonal and 
canonical momentum distributions, (Anfc) CE , versus energy 
per particle, e. (b) Integrated difference of the momentum 
in the microcanonical and canonical ensembles, (Anfc) CE ME , 
versus L~ 119 for e = 1.07J, where a fit to (Anfc) CEME = 
zL -7 gives 7 = 1.0 ± 0.07. Inset: Scaling exponent, 7 of 
(Anfe) CB ME as a function of the energy per particle, e. 



the thermodynamic limit, we calculate the scaling expo- 
nent of the integrated difference of the momentum dis- 
tribution, (An fc ) CE ME = zL-i, where (An fc ) CE ME = 

J2k Iv%)ce - («fc)ME|/Z]fc(^fc)cE- In Fig- 4(b) we plot 
(Anfe) CE ME versus L~~< for energy e = 1.31J, where 
7 = 1.0. In the inset, the exponent 7 is plotted as a 
function of the energy per particle. The mean value is 
7 = 1.03 ± 0.04, which is consistent with L _1 scaling for 
the difference between the canonical and grand canon- 
ical momentum distributions for hard-core bosons in a 
box found in Ref. [l[. 

On general grounds, the canonical and microcanonical 
ensembles are expected to equivalent in the thermody- 
namic limit, which is confirmed by our numerical results. 
We also see increasing agreement between the results of 
the generalized Gibbs and generalized microcanonical en- 
sembles as system sizes increase. There has been signifi- 
cant work in recent years on the equivalence of classical 
ensembles and less in quantum systems. In the classical 
case, Ellis et al. have demonstrated that the microcanon- 
ical and canonical ensembles are equivalent if and only 
if the thermodynamic functions are equivalent, which is 
the case if the microcanonical entropy is concave Q ■ Fur- 
thermore, when the two ensemble are not equivalent, it 
may be possible to construct a generalized canonical en- 
semble, which contains an additional exponential which 
is a continuous function of the Hamiltonian, that is equiv- 
alent to the microcanonical ensemble [3j. 

Scaling of the GGE. Given the trend observed in our 
results, and similar calculations in [J| for much larger 
system sizes, we expect the integrated difference be- 
tween the momentum distribution in the DE and GGE, 
(Arifc)GGE, to vanish in the thermodynamic limit. We 
fitted (Anfc)oGE = zL^ 1 and report the results for 
e = 1.13J in the main text. In Fig. 5 we plot the ex- 
ponent, 7, as a function of the energy per particle of the 
quenched state for all energies studied. The mean value 
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FIG. 5. Exponent of power- law scaling of (Ank)cGE versus 
inverse system size as a function of the energy per particle. 

Number of States. The total number of states in 
our system, the number of states in the ME, and the 
number of states in the GME all scale exponentially 
with the system size as z L . For the total number of 
states, L\j [N\(L — N)\], using Stirling's approximation, 
z w (1 — v) v ^ 1 v^ 1 ' sa 1.65, where v is the filling factor. 
We confirm this scaling by numerically fitting to the to- 
tal number of states. For s = 0.72 J, we find numerically 
that z « 1.51 in the ME, (8 ME = 0.05J) and z w 1.37 in 
the GME (Some — 0.8). Thus the number of states in 
the ME and GME windows as a the fraction of the total 
number of states vanishes in the thermodynamic limit. 
Additionally the ratio of states in GME to ME vanishes. 
These values are typical, although the precise exponents 
depend on the choice of the ME and GME window sizes. 
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FIG. 6. (a) Average density fluctuations, <7 c i C n versus L~°' 5 
for fixed energy per particle in the diagonal (DE), generalized 
microcanonical (GME) and microcanonical (ME) ensembles, 
(b) Momentum fluctuations, o\ for ka — 0. e — 0.72 J, Sme = 
0.05J, 5qmb = 0.8. 
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Fluctuations of Local Observables. For local ob- 
servables, the fluctuation of the eigenstate to eigen- 
state expectation values are expected to scale as L -1 / 2 
[|. In Fig. 6(a) we plot the fluctuations of the 
site occupations averaged over five sites, Cdcn = 
1/5 ELi VW^Im) ~ (hx^N) 2 vs. L-° 5 in the DE, 
GME, and ME distributions along with a linear fit to 
CTdcn versus L~~ 5 . The data strongly suggests that the 
fluctuations scale as L~ 5 as predicted in Ref. [H and 
will vanish in the thermodynamic limit for all three en- 
sembles. 



ensemble it is difficult to reach a conclusion as to whether 
it vanishes or remains finite in the thermodynamic limit. 
We do, however, find that for this observable the results 
are clearly different from those for observables that only 
contain short range correlations, and that the behavior of 
the width of its distribution within the ME behaves quite 
differently from the one within the DE and the GME. In 
the latter two the width decreases much rapidly with sys- 
tem size. 



In the main text, we presented evidence that the fluctu- 
ations of the eigenstate to eigenstate expectation value of 
the momenta occupations vanish in the thermodynamic 
limit for the DE and GME. We have also studied how 
the width of the distribution of the momenta eigenstate 
expectation values, o>, scales for the ME. In Fig. 6(b), 
we plot the fluctuations of fik a =o ■ For the microcanonical 
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